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Discussion of an Uniqueness Theorem 
by L. N. Slobodetsky and M. I. Chramov, 
for Solutions to the Cauchy Problem for a Quasi-Linear 
Symmetrical System of Differential Equations. 

by 

J. K. Billings 



The generalized Cauchy problem is succinctly 
explained by I. G. Petrovsky ([1], page 26). Briefly, 
we are given a system of N equations with N unknown 
functions u^, u 0 , • • • , u^ T 






At 



k k. k 

5x 0 ° 3x, 1 ”* 3x n 






for (i,j = 1,2, ---N) 

and are required to find a solution for u^, u 0 , •••u^ of 
this system in some neighborhood of a smooth n-dimensional 
surface E, this solution to satisfy certain prescribed 
initial conditions. In their paper, ([2]), Slobodetsky 
and Chramov consider the question of uniqueness of a 
solution to this Cauchy problem for particular quasi- 
linear systems of the 1st and 2nd order. Because their 
systems contain a certain type of symmetry, the question 
of uniqueness is answered more easily by a direct 
examination of the system itself, than by an attempt to 
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reduce it to the question of uniqueness of a solution to 
a related linear system, in the manner of Hadamard ([3]). 

We will put most of our attention on the 
symmetric system of the 1st order. We say that the 
system 



^ k 
5 t 



n 



N 



(i) 

u + f k - 1 - 2 >'"> N > 



( 1 ) 



(i) (i) 

(where . = a^,. (t,x,u) and f^ = f k (t,x,u) are complex 
functions of the real variables t and x, and the unknown 
vector — function u jj x = (x^,***, x n ) and u = (u^,***, u^)^ 
is symmetric in some domain D of the (n+N+1) -dimensional 
space (t,x,u), if in D we have 

(i) _(i) 

a jk = (i = = 1,2, “SN) 



( 2 ) 



The domain D, if it contains points M(t,x,u), 
must also contain the points M(t,x,u) where u = (u^ , • • • , Ujj) . 
It is further assumed that the domain is bounded by the 
plane t = 0, and lies in the region where 0 < t, for 
physical considerations. 

Slobodetsky and Chramov showed that a solution to 
system (1) is unique in a certain subdomain of D if the 
following conditions are satisfied: 

1. The system has the symmetry described above. 

2. The functions a are continuous and have 

continuous derivatives in D, with respect to the 
variables x 1 ,**-,x n , u N . 
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3. The vector function f satisfies the 
inequality 

|f(t,x,u) - f(t,x,v)| < 4(|u-v|) (3) 

when (t,x,u) and (t,x,v) ^ D. 



4. F(Z) = 1/ Z 8i(Z) is a positive, increasing, 
and concave function, such that 

S 



f diverges, 0 < $ 



(4) 



The quantity |u-v| is defined as follows: 

|u-v| = (u-v, u-v) 1/0j 

where the scalar product of two vectors Z = (Z^,* #, ,Zj T ) 
and y = (y^***^) is given by 

N 

(z,y) = 2] z i y-j* 






Thus 



u-v =. 




1/2 



3= 1 



(u r v o> vV' 



F(Z) is called a concave function in some 
interval (a,b) if for any Z k £(a,b) (k = 1, 2,***,m) 
and for any non-negative 1^ (k = 1, 2, such that 

m _ 

1,^ = 1, we have the inequality 



k=l 

m 



\ 



k=l 



k s H kk 



k=l 



/ 



It is remarked that these conditions of 
uniqueness of the solution for such a system appear 
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similar to the conditions of Osgood ([4], page 344), 
as discussed by J. Tamarkine ([5]) for the uniqueness 
of the solution of the Cauchy problem for a system of 
ordinary differential equations. The condition of 
Lipschitz for the uniqueness of a solution to the 
equation = f(x,y) is the well known inequality 



In his paper ([4]), Osgood extended the Lipschitz 
condition in the following manner. Given a continuous 
function cj (u) such that 



(b-y^l < S, |b-y 0 | < S, then at the point x = a, the 
solution will take the value b and this solution will be 
unioue. Tamarkine reformulated this condition by noting 
that if 



where it) is a continuous, positive, increasing function, 
and u) (0) = 0, then the solution is unique if, 



|f(x,y 2 ) - fCx^)! < k I y 2 -y x I • 



CO( 0) =0; 0 < £P(u), 0 < u; <V( -u) = o^(u) 



u 

and lim / 

£ -»o C 




then if 



|f(x>y 2 ) “ f I ^ w(y 2 -y 1 ) for |a-x| <f> 



|f(x,y 2 ) - f^ypi < c^dyg-yj^l) 



for 0 < u < u , 




We shall now attempt to demonstrate the proof 



that these conditions stated by Slobodetzky and Chramov 



do indeed insure uniqueness in the very particular case 
when n = 1 and N = 2. The system (1) can then be 
written out explicitly as 
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^ U^ 




d u 0 




at - a n 


ax + a 2i 


sir + 


f x (x,t,u) 




c* u JL 


^U 9 




at “ a i2 


Sx + a 22 


ax + 


f ^(x, t , u) 



(5) 



which of course may be also written in matrix form 



(M- (a) (t^) + < f > 



where A 




x2 




( 6 ) 



This symmetrical system (5) is of the hyperbolic type 
described by Petrovsky ([lj page 53) because its 
characteristic equation 

- >5^1 = 0 (7) 

has only real roots for any real <=»< . This is readily 
seen by expanding the determinant to get 



V - (a 
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+ a 22^^ “ ^ a ll a 22 + a 12 a 21^° < ' = 0 



By the symmetry property assigned to this system we note 

q 

that a^ and a^ are real, and that a-^ a 0 -^ = la^l • 

2 2 2 

Thus the discriminant is e< (a^ra 0 >) + Ha^^l 3 

which, for a real , is always greater than zero. 

Since no restrictions have been placed on the multiplicity 
of roots of the characteristic equation (7), this system 
does not, generally speaking, belong to the class 



"hyperbolic in the narrow sense" discussed by 
Petrovsky ([1] page 59). 
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Because of the symmetry property assigned to (6), 
the matrix A is Kermitian, or self-adjoint, and there is 
associated v/ith it the Ilermitian form 

H(y:<*) = (A<s<y,y) (8) 

N 

V 

where (A<=<y,y) = s' (A^y) . y. is again the scalar 

j=l 3 3 

product of the vectors A«<y, and y. Since the a^ are 
to be continuous on a closed domain to be described below, 
we can say that I a jj c l < <1 a nd I j^l < for = 1, where 
^ are the roots of equation (7). It is easily shown 
that L < Mq, where M is an absolute constant independent 
of the coefficients of system (6). We may further point 
out a property of the Ilermitian form, 

|H(y: cS-') ^ L \y\ 2 ^< (9) 

where |y| = V (y^y) is the length of the vector y. 

The domain _n_(t) in the 2-dimensional space 
(t,x) is called "fundamental" for equation (6) if it is 
restricted to be between t = 0 and 0 < t = T and is 
bounded by a smooth curve E, the exterior normal to 
which forms an acute angle (n,t) with the t-axis, such 
that tan(n,t) < . We further denote by S(t) any 

section of_n_(T) by a line t = constant. See figure 1* 

V/e will show that for any solutions u(t,x) and v(t,x) 
which satisfy the prescribed initial conditions, and 
such that (t,x) £ _n_ (T) => (t,x,u) and (t,x,v) £ D, 



+ 




Figure •//!• 



bicetch of X1(T) 



then u(t,x) £ v(t ,x) in the fundamental domain .fl(T). 

Since u and v are both assumed to be solutions 
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of (6), we may subtract one from the other to obtain 



or 



= A ( t > x ’ u )|“£ + f(t,x,u) - A(t,x,v) ^ - 
^A(t,x,u) - A(t,x,u)|^J 

^ = A(u)^l + [ A (u)-A(v)] |X + f(u) -f(v) 



a x 

- f(t,x,v) + 



( 10 ) 



where for brevity we write A(t,x,u) = A(u), A(t,x,v) = 

A(v), f(t,x,u) = f(u), and f(t,x,v) = f(v). Multiplying 
equation (10) from the right and from the left by the 
vector (u-v) and adding the two results we obtain 

, . (ID 

£| U -vl 2 

A(u)-A(v)J|j ) o-v^ + 2 Re^f(u)-f(v),u-v) 

The last two terms on the right hand side of (11) are a 
result of the Hermitian type of symmetry assigned to the 
matrix A, and the definition of the scalar product employed. 

If we now take the total derivative of the 
expression ^A(u)(u-v), u-vjwe see that 



+ 2 re 



d 



*\ 



(uu)(u-v), u-v)|' = (u-v), u-v)+ 

+ ^A(u) u-v)+ 

+ ^A(u) (u-v) , 



( 12 ) 
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But ^A(u)(u-v), - = 



'a^u) 






a 12 (u > 




- (U 1' V 1 ) 



+ (u 2 -v 2 ) 



M u l" v l-> CIV-2 '2 

a ll (U > + a "’ ( “ ) 



, a(u 2 -v,) 

12 '“) b 



a 21 (“) 



3 < u l- v l> 



- a 22 (u > 



a(u 2 -^7 



while 



f (u-v), A(u) ^ . O L ~ z! ) = 
\ d x 7 



= f(“l- v l U 2 -V 2 



(u l' V l) 



! (U 2 -V 2 ) 




— ^ a(u 1 -' , 1 ) , „ mu 2 -v 2 ) 

a 12 ?>x + a 22 ( “) ■Jx 



so that (a(u)(u-v), ^ 4 U_V > A ( u ) / 

and hence (12) may be rewritten 

( A(u) A^vl, u -v) + (u-v, A( u) ISSXl). 

(a(u)(u-v), u-v - d j^“) (u-v), u-v) 



( 13 ) 



d 

dx 
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From the condition of continuity imposed on the 

d A 

derivatives of a.. , we know that the matrix is 

jk dx 

bounded in the fundamental domain so that we may assert 
that 



(u-v), (u-v>) | < |u-v| 2 where is a 



(14) 



positive constant. 

Now by expanding the matrix product and noting 
that the absolute value of a sum (or difference) is less 
than or equal the sum of the absolute values of the 
individual terms, we see that 



J Re ^(A(u) - A(v)j -|-£, u- 

^ j a ll (u > -53E + Re a 21 (u > Tx 

I d v i V 

+ ) a ii (v > TZ + Re a 2i< v > S 



u i- v l' 



1 b v i b v 2 

+ j Re a 12 (u) + a 22 (u) 

+ | Re a 12 (v) + a 22 (v) |-§ 



u l“ V l' + 



|u 2 - v 2 I + 



|u 2 -v 2 I . 



Again utilizing the boundedness of the elements a^ c and 
the derivatives involved, we see that 

b v 



2| Re i^A(u) - A(v)^ 



b x 



, u-v 



( 15 ) 



b |ui- v i| + C| u 2-v 2 I <M 2 |u-v| 2 
where B, C, and M 2 are positive constants, 
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From (3) we note that 

2 1 Re ^f(u) - f(v), u-vy | < 2|u-v| $ (|u-v|) = 

= 2f(|u-v| / 2 . (16) 

Then substituting (14), (15), and (16) into (11) we 
see that 



d I u-v 1 



< 



a t s dx 



^A(u) (u-v) , u-v) + 



+ | u-v | 2 + 2F ^ | u-v 1^ 



(17) 



where is a positive constant. Let us now integrate 
(17) over the fundamental domain _n_(t) for 0 < t < T, 
by applying Green's formula. Taking this term by term, 
and observing that u H v for t = 0, we obtain the 
following: 



2 

/ -^ - y V -l — dt dx = / | u-v | 2 dx = 

-TU ^ C 



= / t u-v | 2 dx + / t u-v | 2 n , dE 
S(t) E 



where n^ = cos (n,t) . 



• ^A(u) (u-v) , u-v) 



dt dx = 



= / asr? H(u-v:l)\dt dx = / H(u-v:l)n v dt = 



= f H(u-v:l)n dt + / H(u-v:l)n dt 

X ^ A 

S(t) E 

where = n = sin(n,t) = cos(n,x). 

A. 



= / H(u-v:®<)dt. 

E 
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• / |u-vr dx + / | 

• • S(t) E L 



- H ^u-v: cos(n,x)J 



| u-v| ^ cos(n, t) - 



(18) 



dE < M_ / | u-v| 2 

= 3 0 S(r) 



+ 2 / d *3*“/ Fdu-vl^) dx. 

0 S(» 

According to inequality (9) we may assert that 
|u-v| 2 cos(n,t) - H(u-v:cos(n,x)) > |u-v| 2 |^cos(n,t) - 
-L sin^n > t)J > 0, since L < ctn(n,t) = 

Thus we have 

/ |u-v| 2 dx < M, f d l' f |u-v| 2 dx + 



S(t) 



0 



SO) 



+ 2 / d'J* f F(|u-vr) dx. 

o s(r) 

Because F(Z) is concave, we know that in an interval 



(19) 



(a,b), where Zj,€.(a,b) and 



X k " 1 



m 



m 



g l k F(2 k ) < rg X k Z^ 

Multiplying and dividing by the length of the interval, 
say b - a, gives 

■RT 2 l k P(Z k )(b-a) S F Lis- )X W-“) 

k=l \ (Ki.l fj 

Then if we let m increase indefinitely, we see that 

l k (b-a) becomes a k-th subinterval of (b-a) and each 

side is a Riemann sum, so that 



w / F < z > ax s f (-eV i H- 



( 20 ) 
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Applying (20) to equation (19) we see that 

/ |u-v| 2 dx < M- / d O'/ | u-v| 2 dx + 

s(t)_ = 3 o s((H 

+ 2 / o- Wf /~ / |u-v| 2 dx j d^- (21) 

0 V 3 ^' S(0“) / 

where cr(7) is the length or measure of the line segment 

sm. 

Without loss of generality we may assume that 
< 5 ~ (fY) <1. Again taking into account the fact that 
F(Z) is concave, we obtain that 



o-(0')F(x 1 ) + (l-<r(0-)) F(x 2 ) ^ F ( x 1 cr( T)J + x 2 (W(?“)y . 

If now we let x^ = 0, and observe that F(0) = 0 we note 
that 

<5-C»F(Z) ^ f(cp(t)Z.) (22) 

If now we let Z(t) = / |u-v| 2 dx and apply (22) to 

S(t) 

equation (21) we find 



Z(t) ^ M 3 / Z(7-)d y + 2 / F(zm) d'T- (23) 

and the proposition will be proved if we show Z(t) 5. 0 
for 0 < t < T. Let us suppose the contrary. We know 
that Z(0) = 0. Therefore there must exist two numbers 
o< > 0 and ft > 0, such that 



(a) Z(t ) = 0 when 0 < t < . 

(b) 0 < Z(t) < 1 when c< < t < ^. 

Then for c< < < /3 , we see from equation (22) that 

F(l)Z( , 3 k ) < F/Z(T)) , and hence 
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Z(t) ^ 




d'T + 



+ 2 




dr 




dr, 



or Z(t) ^ M 4 / f(z (r)J dr , o< < t < P , (24) 

where since F is an increasing function, F(l) > 0. 



Now put § (t) = / F^ZCy)J d?- . 

an increasing function we have: 



Then, since F(Z) is 



= H 4 F(z(t>) < m 4 F (§<t) ) . 

h?) d < 

Evidently / 

which contradicts the divergence of the integral 

s 



f 



dS 



0 ^ 



, 0 < £ . QED. 



(25) 

(26) 



Thus the uniqueness has been demonstrated. 

To demonstrate the practical applicability of 
these uniqueness conditions let us consider briefly 
Maxwell's wave equations. Of course, solutions to 
Maxwell's equations have been attained in other ways, 
(c. f. Courant and Hilbert [6]), but this will serve to 
illustrate our principles. Maxwell's wave equations 



are : 
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(1) V X H = i 

(2) V x E = - i -|| 

(3) V . H = V .E = 0 

where E = E(x,y,Z,t) is the electric field strength 
vector, H = H(x,y,Z,t) is the magnetic field strength 
vector and C = constant velocity of light in vacuum = 

3 x 10^ cra/sec. 

Obviously for our purposes we need consider only 
system (1). In component form we have 

BE, oH, 3 H 0 

- 1 = C C — 

B t ^ x 2 ^ x 3 




BH, 

0 IT 1 - " c T - 2 

O X3 0 



o E- BH 0 ^H, 

±- = C C — . 

d t B x^ ^ x 2 

Denoting E^ , E ? , and E^ by u^, u^, and u^ respectively 
and H^, Hg* by u^, u^, Ug respectively, Maxwell's 
equations take the form of system (1). We see that our 
coefficient matrices will then be as follows: 
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( 1 ) 



(1) 



A = (a jk ) 



( 2 ) 



( 2 ) 



A = (a jk ) 



(3) 



(3) 



A = (a jk ) 



and in matrix form, 




0 0 \ 



0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 -C 

0 0 0 0 0 0 

0 0 0 0 0 0^ 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 c 

0 0 0 0 0 

0 0 0 -C 0 0 

0 0 0 0 

0 0 0 0 0 




i=l - - 1 

If now we multiply both sides by V -1 I our matrix A 
takes on the desired Hermitian type of symmetry and we 
are ready to test our conditions. 

1. The system does have the required symmetry. 

2. The functions a , being constants, are 
continuous and do have continuous derivatives in any 
prescribed domain, with respect to all variables. 



(27) 
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3. Since f, c 2 0, the inequality 

|f(t,x,u) - f(t,x,v)| < $ (|u-v|) 



is satisfied for $ (|u-v() = |u-v|, say. 



4. Then F(Z) = fz $ Cf~Z~) = Z is a positive 
increasing concave function such that 

S 

T(zT 



r dZ r 

J J 



dZ 



0 



(0 < S). 



Although in the case of Maxwell’s equations as 
stated, the vector-function f is identically zero, we 
note that the uniqueness of the solution is not 
destroyed even though a perturbation effect is present, 
so long as the perturbing function satisfies the 
inequality conditions of 3 and 4. 

We will conclude with a statement of the con- 
ditions required for uniqueness of the solution of a 
symmetrical system of second order, which in matrix 
notation has the form 



a 2 u 

3t 2 




^ 2 u 

3 x i axj 



+ 



+ 



f(t,x,u, 






(28) 



The proof of this uniqueness can be found in reference 
work [2j, and is very similar to the proof for the first 
order system. 
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The symmetry property appears in this statement, 
that the matrix 




j (t ,x,u, 



^u 

-&t ’ 



3 u 
B x. 



~B u 
D x 



n 



) 



is Hermitian, and A i ^ = A^. It is further assumed 
that for any n-dimensional vector ^ , 




(29) 



where 0 is a real number. From (29) follows the 
fact that for any real <=< the characterisiic equation 



n 




i, 3=1 



A. . 

i 3 




0 



has only real roots, and hence system (28) is hyperbolic. 
Let the following conditions be satisfied: 

1. The matrix A is continuous and is 
continuously differentiable with respect to the variables 



x 



1’ 



9ui 

W‘”’V“aT’ 



. i^N j^i 

’ -5t ’ 3x x ’ 



3 u 



N 



^ X n 



2. The vector function f satisfies the 

inequality 



f(t,x,u, 



3 u 3 u 
3t’ dx x ’ 



d u 

3x n 



) - 



- f (t , x, v, 



~3 v 

3t’ 





< 
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< S 4 



where 







n 


u- 


■v| 2 + 


z 






k=0 




u 


~d u 




X 0 " 


d t 



'b x, 



k 




3. F(Z) = V Z $ (ITT - ) is a positive, 
non -decreasing, and concave function such that 

s az 

/ p - ^ - gy diverges (O < S )• 



Then if u and v are two twice-continuously 
differentiable solutions of system (28) which satisfy 
the initial conditions 



u 



t=0 



t=0 



~b u 
c* t 



t=0 



3 V 

at 



t=0 



these solutions must coincide in some fundamental domain 
of system (28). As a final note, this uniqueness theorm 
is valid in particular, for a single quasi-linear 
hyperbolic equation of the second order. Thus the work 
of S. L. Soboleff [7] in regard to obtaining a priori 
values for single hyperbolic equations of second order 
can perhaps be extended analogously to such symmetrical 
systems. 
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